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Abstract We present a method for a complete char- 
acterization of a femtosecond ultraviolet pulse when a 
fundamental near-infrared beam is also available. Our 
approach relics on generation of second harmonic from 
the pre-characterized fundamental, which serves as a ref- 
erence against which an unknown pulse is measured us- 
ing spectral interference (SI). The characterization appa- 
ratus is a modified second harmonic frequency resolved 
optical gating setup which additionally allows for tak- 
ing SI spectrum. The presented method is linear in the 
unknown field, simple and sensitive. We checked its ac- 
curacy using test pulses generated in a thick nonlinear 
crystal, demonstrating the ability to measure the phase 
in a broad spectral range, down to 0.1% peak spectral 
intensity as well as retrieving 7r leaps in the spectral 
phase. 
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1 Introduction 



Techniques allowing for a complete characterization of 
the femtosecond pulse instantaneous intensity and fre- 
quency were first introduced in the nineties. The FROG 
(Frequency Resolved Optical Gating) [1,, SPIDER (Spec- 
tral Interferometry for Direct E-field Reconstruction) [2] 
and sonographic techniques [3] are the most popular. 
These techniques are usually based on sum frequency 
generation of an unknown near infrared (NIR) pulse with 
a modulated copy of itself. The resultant UV radiation is 
registered as a function of modulation. Finally, the pulse 
envelope and phase are retrieved from acquired data. It 
is crucial for those methods to filter out the fundamental 
components from the frequency sum signal. Fortunately 
it can be easily accomplished by inserting a color glass 



* tel: +48 56 6113332, fax: +48 56 6225397, E- 
mail: mkacpr@gmail.com 



Fig. 1 Pictorial diagram of the measurement scheme. £ U nk'- 
unknown pulse, £f: fundamental pulse, £ Ic f. reference second 
harmonic, X: thin nonlinear crystal. 



filter. Moreover the UV signal is measurable with inex- 
pensive silicon-based photodetectors. 

Frequently one faces the need to characterize the sec- 
ond harmonic (SH) of the pulse of ThSapphire lasers. 
This is the case in the ultrafast spectroscopy [HE], mi- 
cromachining [5] or downconversion-based photon pair 
sources [TUgpiTO] . Directly generalized FROG and SPI- 
DER techniques suitable for characterization of the UV 
pulses are based on upconverting or downconverting the 
unknown second-harmonic pulse using the fundamental 
beam |llU12j . The resultant radiation, which carries in- 
formation about the unknown pulse, is difficult to de- 
tect. In the case of using upconversion it is centered 
around 266 nm where the silicon detectors are ineffi- 
cient. On the other hand using downconversion produces 
weak pulses at the fundamental wavelength that are eas- 
ily overwhelmed by the stray background. Moreover the 
nonlinear up- or downconversion process requires sub- 
stantial intensity of the measured beam. The latter re- 
quirement has recently been diminished by application 
of an optical parametric amplification based conversion 
process [T51ll4| , however these new approaches are com- 
plicated in use. 

In this paper we present a method for complete char- 
acterization of SH pulses when the fundamental beam is 
also available. This is a typical experimental situation, 
which can be exploited to avoid nonlinear conversion of a 
SH beam. It is assumed, that the fundamental pulses are 
characterized prior to the measurement, for example us- 
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ing the well-established second harmonic generation fre- 
quency optical gating (SHG FROG) technique [1 . The 
scheme of our method is presented in Fig.[T] A portion of 
the fundamental beam is converted into reference second 
harmonic. Thanks to application of a very thin crystal in 
this step the envelope of the reference pulse can be pre- 
cisely calculated. Then the reference and the unknown 
pulse are brought to interfere on a slit of a spectrom- 
eter and spectral fringes are registered. The phase of 
the unknown pulse is retrieved from this signal using 
the Fourier filtering technique |15j . We demonstrate this 
method for a pulse generated in 1 mm thick beta barium 
borate (BBO) crystal oriented for type I phase match- 
ing. We were able to reconstruct the phase in spectral 
range where the pulse intensity is above 0.1% maximum, 
including the it phase leaps. 

2 Method 

The first step in our measurement method is registering 
FROG trace for the fundamental beam pQ. Next we run 
the standard retrieval procedure using FROG software 
Femtosoft Technologies) and obtain a complete informa- 
tion on the electric field of the fundamental pulses. In 
particular we learn about their spectral envelope £f(oj)' 



£ f (lo) = v^Re'^M, 



(1) 



where If{oj) is the spectral intensity and 4>f(w) is the 
spectral phase. Note that Jf(w) can be measured directly 
with a spectrometer. This gives us a possibility to verify 
the accuracy of the data retrieved from FROG. Next 
we double known fundamental pulses in a BBO crystal. 
The crystal is thin as compared to the distance on which 
NIR and UV pulse diverge in time by their duration due 
to crystal dispersion. We take advantage of this fact to 
calculate precisely the output second harmonic and use 
it as a reference. It's spectral envelope £ le f(ui) is given 
by the usual convolution formula: 



£ re f(to) cx / duo' £f{oj')£f{oj — w'). 



(2) 



In particular we learn about the spectral phase of the 
reference pulses (f> rc {(oj). Note that above formula is a 
reflection of the fact that the second harmonic field in 
time is the fundamental field squared: £ re f(t) — £f{^)- 

We exploit the spectral interference method in order 
to retrieve the difference of phase between the unknown 
second harmonic pulse £ U nk(w) and the reference pulse 
£ r ef(ct>). This is accomplished by directing them with a 
relative delay r into the spectrometer slit, where they 
interfere. Hence we measure the interference spectrum 
Isi(oj) of the form: 

I — ' 1 2 

Isi(u) = \£ Te f (oj) e~' lUJT + f UI1 k M | 

= J re f (oj) + 7 unk (<j) 



-2y/l Ie { (w) J unk (oj) 
X COS (<t>rei(u) ~ <?WM - Wt) . (3) 



An exemplary interference spectrum Jsi(w) is shown in 
Fig. [2j The key information about the unknown phase 
</>unk(w) is contained in the cos (4> le {(u>) — 0unk(<J-O — ojt) 
term. First, we retrieve </> re f(<x>)— </> U nk(w) using the Fourier 
filtering technique [15 . The interference spectrum can be 
separated into three distinct elements: the sum of spec- 
tral intensities J ro f (w) + J un k {oj) which is slowly varying 
with oj and the fringes described by the last term in |3| , 
which consist of positive and negative frequency parts 



varying as e 



and e 



We calculate the Fourier 



transform of the interference spectrum, which splits those 
parts in the Fourier domain, provided the delay r is 
large enough. A typical result is shown in Fig. [3j Next, 
the Fourier transform is multiplied by a supergauss filter 
function F(t): 



F(t) = exp 



At 



(4) 



where t denotes time in the Fourier domain and At is 
the filter width. F{t) is nonzero near the t — t point and 
zero elsewhere. Thus the product of the filter function 
F(t) and the interference spectrum in the Fourier do- 
main Isiyt) contains only the term originating from the 
part of Jsi(w) which varies as e +luJT . It is transformed 
back into the frequency domain, yielding: 

4l M = V ^rrf M -?unk (w) 

X exp [i OrefM - </>unkM + Wr)] (5) 

The argument of Jgj (u>) is the phase difference U nk(^)~ 
</> re f(w) and the linear term ojt. Hence the unknown phase 
is given by the following formula: 



v unk 



arg 



t) Ic{ + LOT. 



(6) 



Finally the spectral intensity of the unknown pulse J un k(w) 
is measured with a spectrometer and the spectral enve- 
lope can be reconstructed, which concludes the charac- 
terization: 



funkM = V^unkHe l0U " kM - 



(7) 



Note, that in the experiment nearly identical pattern 
of the fringes is obtained for a given magnitude of the 
relative delay r between the pulses, regardless of its sign. 
On the contrary, the above formulas remain valid only 
when a correct sign of r is used: from the convention 
used in Eq. (|3| it follows that if the reference beam path 
is delayed then positive r is to be inserted. Also, note 
that a small variation in the assumed value of r, diffi- 
cult to avoid in any real experiment, results only in an 
irrelevant temporal shift of the reconstructed unknown 
pulse shape. 
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Fig. 2 Interference fringes measured by the spectrometer 
used in the method. For a detailed description of the setup 
see Sec. [3] Arrows mark points of ir spectral phase leaps. 
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Fig. 3 The Fourier transform of the interference fringes pat- 
tern (thin line) and the filtered signal (thick) obtained using 
the supergauss filter form Eq. [4] with filter width At — 500 fs 
and temporal delay r = 1000 fs. Note the logarithmic scale 
on vertical axis. 
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Fig. 4 Experimental setup. BS1, BS2: beamsplitters; FM, 
CCM: mirrors; LI, L2: lenses; X: SHG crystal; D: diaphragm; 
BF: blue filter. 



FM. Both thus produced pulses are focused on a 0.05 mm 
thick BBO type I crystal X where the second harmonic 
generation occurs and three beams emerge. Two of them: 
upper and lower are blocked by a diaphragm D. The 
middle beam is filtered out of a scattered fundamental 
beam by a blue filter BF and focused with a lens L2 
on the spectrometer slit. We retrieved the fundamental 
pulse with a FROG error [T] of 0.0012 and we checked 
that the retrieved spectrum is consistent with the fun- 
damental beam spectrum measured directly by the spec- 
trometer. 

The second step is to measure the interference of the 
reference pulse and the unknown UV pulse generated 
by an experiment. For this purpose we switch the posi- 
tion of the flipping mirror to direct the generated UV 
pulse to the spectrometer. This pulse is assumed to pass 
unaffected through the crystal X, while the NIR beam 
reflected off the cornercube mirror CCM generates sec- 
ond harmonic of the form given by Eq. [2] These two UV 
pulses are let through by removing the diaphragm D and 
made interfere on the slit of the spectrometer. 

Note that the crystal suitable for FROG must be thin 
enough to guarantee the validity of Eq. [2] Therefore we 
can use the same crystal for reference pulse generation 
as for SHG FROG measurement. 



Our measuring apparatus is depicted in Fig. [4] A 780 nm 
centered beam from a mode-locked Ti:sapphire laser is 
split on a beamsplitter BS1 into two paths. One part 
goes directly to the FROG-SI measurement apparatus, 
while the other is sent to the test pulse preparation arm. 

In the first step we measure the SHG FROG trace of 
the fundamental beam. This measurement requires two 
replicas of the fundamental beam, which are produced 
by the beamsplitter BS2. One replica reflected off BS2 
undergoes an adjustable delay by bouncing off the cor- 
nercube mirror CCM, while the second one, transmitted 
through BS2, is directed with help of the flipping mirror 



4 Results 

We verified the accuracy of our method using test pulses 
generated in a 1 mm thick BBO crystal and filtered out 
of remaining NIR light and attenuated with neutral den- 
sity filters. The crystal and filters are represented as an 
"experiment" in Fig. [4] The test pulse is directed to the 
FROG-SI apparatus. When producing the test pulses 
we did not focus the fundamental beam in the crystal to 
achieve a well defined phase matching angle. This simpli- 
fied the calculation of the test pulse envelope and made 
the spectral phase features more pronounced. 
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In Fig. [5] we plot a comparison of a theoretically 
calculated pulse and one retrieved experimentally. The 
results displayed in Fig. [5ja) were obtained for thick 
crystal oriented for phase matching at the wavelength 
778 nm, which is equal to the central wavelength of 
the incident fundamental beam. Fig. ^b) shows the re- 
sult of a similar measurement with the crystal tilted by 
63 mrad, which results in phase matching at 785 nm. 
Note the fidelity of it phase leaps at wavelengths where 
the spectral intensity meets zero. The leaps can be no- 
ticed in raw interference fringes Fig. [2] The inversion 
procedure naturally fails at wavelengths where the un- 
known pulse intensity is below the noise level and the 
fringes disappear. From Fig. [5] we infer that the thresh- 
old spectral intensity above which phase reconstruction 
is possible is about 0.1% of the maximum. Note that 
our method is capable of retrieving phase in case when 
the spectrum of the unknown pulse is nonzero in two or 
more separated regions. 

We calculated the pulse intensity and phase at the 
output face of the crystal producing test pulses using 
the formula for second harmonic £sn(u) generated in a 
thick nonlinear medium in case when pulse distortion 
can be neglected [IB] : 



£ S h(w) oc sine - 2uj pm ) 



x I duj' Ai(uj')A 2 {uj - ui') 



(8) 



where L is the crystal thickness, u> pm is the frequency 
of perfect phase matching and A/3 is the difference of 
the reciprocals of group velocities of fundamental and 
second harmonic waves. 

The necessary presence of neutral and color glass fil- 
ters in the test pulse path influences it by contributing 
an additional phase A<j>&it ei (iv) . We describe it by Taylor 
expansion around the pulse central frequency loq: 



A(j) mteI (oj) = f3 f+f3if(uj- 



(9) 



The constant and linear phase terms (3q f and Pi / corre- 
spond to pulse retardation. The only relevant and appre- 
ciable contribution is the quadratic dispersion /3 2 //2(w — 

For measuring 02 f we replaced a 1 mm BBO crys- 
tal with one that was 0.1 mm thick, which generated 
more broadband test pulses. Next, we acquired a refer- 
ence spectral interferogram. Then we inserted additional 
filters, identical to those to be characterized, in the test 
path and registered the interference again. The difference 
of the spectral phases retrieved in those two cases equals 
A(f> futeriy) ■ By fitting the latter by a square polynomial 
we computed (3 2 f- Using these data we could substract 
the influence of the filters from the measured test pulse 
phase obtaining the result shown in Fig. [5] 

The agreement of measurement results with theoreti- 
cal models certifies the validity of the presented method. 
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Fig. 5 Spectral amplitude (solid line-simulated, rectangles 
- retrieved experimentally) and phase (dashed line - simu- 
lated, circles - retrieved) obtained for thick crystal oriented 
for phase matching at 778 nm (a) and with the crystal tilted 
by 63 mrad which results in phase matching at 785nm (b). 



5 Summary 

We have demonstrated a novel method for a complete 
characterization of the near-UV second harmonic of ul- 
trashort pulses. It is applicable in situations when an in- 
tense NIR beam is available. Additionally, the unknown 
pulse spectrum must lie within the spectral range of sec- 
ond harmonic generated in a very thin crystal. 

The main advantage of the presented method is a 
simple setup. It is a hybrid of standard SHG FROG 
with a spectral interferometer and there is no need to 
employ additional spectrometer or crystals. The pulse 
shape is inferred from an UV spectrum which is easily 
registered. The algorithm used for phase retrieval is sim- 
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pie and non-iterative. In addition, our method is linear 
in the unknown pulse field. 

The interference acquisition is single shot, without a 
need to scan a range of time delays r. Additional aver- 
aging over various r is possible with compensation of a 
linear phase. 

Finally, our method allowed us to characterize rather 
small phase contributed to the UV pulses by glass filters. 
Measurement of such a small phase would be very diffi- 
cult to achieve using FROG. 
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